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ABSTRACT 

Context. The Concordia base in Dome C, Antarctica, is an extremely promising site for photometric astronomy due to 
the 3-month long night during the Antarctic winter, favorable weather conditions, and low scintillation. 
Aims. The ASTEP project (Antarctic Search for Transiting ExoPlanets) is a pilot project to discover transiting planets 
and understand the limits of visible photometry from the Concordia site. 

Methods. ASTEP South is the first phase of the ASTEP project. The instrument is a fixed 10 cm refractor with a 4kx4k 
CCD camera in a thermalized box, pointing continuously a 3.88 x 3.88°^ field of view centered on the celestial South 
pole. We describe the project and report results of a preliminary data analysis. 

Results. ASTEP South became fully functional in June 2008 and obtained 1592 hours of data during the 2008 Antarctic 
winter. The data are of good quality but the analysis has to account for changes in the PSF (Point Spread Function) 
due to rapid ground seeing variations and instrumental effects. The pointing direction is stable within 10 arcseconds on 
a daily timescale and drifts by only 34 arcseconds in 50 days. A truly continuous photometry of bright stars is possible 
in June (the noon sky background peaks at a magnitude R ~ 15arcsec~^ on June 22), but becomes challenging in 
July (the noon sky background magnitude is R 12.5arcsec~^ on July 20). The weather conditions are estimated 
from the number of stars detected in the field. For the 2008 winter, the statistics are between 56.3 % and 68.4 % of 
excellent weather, 17.9 % to 30 % of veiled weather (when the probable presence of thin clouds implies a lower number 
of detected stars) and 13.7 % of bad weather. Using these results in a probabilistic analysis of transit detection, we 
show that the detection efficiency of transiting exoplanets in one given field is improved at Dome C compared to a 
temperate site such as La Silla. For example we estimate that a year-long campaign of 10 cm refractor could reach an 
efficiency of 69 % at Dome C versus 45 % at La Silla for detecting 2-day period giant planets around target stars from 
magnitude 10 to 15. The detection efficiency decreases for planets with longer orbital periods, but in relative sense it 
is even more favorable to Dome C. 

Conclusions. This shows the high potential of Dome C for photometry and future planet discoveries. 
Key words. Methods: observational, data analysis - Site testing - Techniques: photometric 



1. Introduction 

Dome C offers exceptional conditions for astronomy thanks 
to a 3-month continuous night during the Antarctic win- 
ter and a very dry atmosphere. Dome C is located at 
75°06'S - 123°21'E at an altitude of 3233 meters on a sum- 
mit of the high Antarctic plateau, 1100 km away from the 
coast. After a pioneering summer expedition in 1995, the 
site testing for astronomy begun in the early 2000's. It 
revealed a very clear sky, an exceptional seeing and very 
low wind-speeds (|Aristidi et al.|[2QQ3l . l2QQ5l : iLawrence et alJ 



120041 : lAshlev et al.ll2005bl : iGeissler Masciadril[2006h . The 
French- Italian base Concordia was constructed at Dome C 
from 1999 to 2005 to hold various science experiments. 
Sum mer time astrononi y experiments have been carried out 
(e.g. iGuerri et al.ll2QQ7l ). The study of Dome C for astron- 
omy during night-time has considerably expanded since the 
first winter-over at Concordia in 2005. The winter site test- 
ing h as shown an excellent seeing abov e a thin b oundary 
layer (|Agabi et al.ll2006l : iTrinquet et"al] [2008: Aristi di et all 
I2009D . a very low scintillation fKenvon e t al.ii2006 l and a 
high duty cycle (Mosser fc Aristidi. 2007.) . Low sky bright- 
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ness and extinction are also expected ([Kenvon fc StorevI 
120061) . 

Time-series observations such as those implied by the 
detection of transiting exoplanets should benefit from these 
atmospherical conditions and the good phase coverage. This 
could potentially greatly improve the photo metric precision 
when compared to other temperate sites (Pont fc Bouchv 



2005 ). A first photometric instrument, PAIX (Chadid et al. 
20071 ). was installed at Concordia in December 2006. A 



lightcurve of the RR Lyrae varia ble star Sara o v er 16 nights 
in August 2007 is presented in IChadid et al.l (|2008f ). and 
results of the whole campaign from June to August 2007 
have been submitted. The sIRAIT instrument also obtained 
light curves over 10 days on the stars V841 Cen an d VI 034 
Cen ([Brigudio et al.ll2009l : IStrassm eier eralll2008 h. 

The ASTEP project (Antarctic Search for Transiting 
ExoPlanets) aims at determining the quality of Dome C as 
a site for futur e photometric surv eys and to detect tran- 
siting planets ([Fressin et al.l l2005l ) . The main instrument 
is a 40 cm Newton telescope entirely designed and built 
to perform high precision photometry from Dome C. The 
observations will start in winter 2010. A first instrument al- 
ready on site, ASTEP South, has observed during the 2008 
and 2009 winters. 

We present here the ASTEP South project and results 
from the preliminary analysis of the 2008 campaign. We 
first describe the instrument, the observation strategy and 
the field of view. Section [3] discusses the main features ob- 
tained when running this simple instrument from Dome C: 
influence of the Sun and the Moon, PSF and pointing vari- 
ations, as well as temperature effects. In section S] we detail 
our duty cycle and infer the weather statistics at Dome C 
for the 2008 winter. These results are combined to a prob- 
abilistic analysis to infer the potential of ASTEP South for 
planet detection and to evaluate Dome C as a site for future 
planet discoveries. 



2. Instrumental setup 

2.1. The instrument 

ASTEP South consists of a 10 cm refractor, a front- 
illuminated 4096x4096 pixels CCD camera, and a simple 
mount in a thermalized enclosure. The refractor is a com- 
mercial TeleVue NPlOl and the camera is a ProLines se- 
ries by Finger Lake Instrumentation equipped with a KAF- 
16801E CCD by Ko dak. For the choice of the camera see 
ICrouzet et al.l ([2007 ') . Its quantum efficiency peaks at 63 % 
at 660 nm and is above 50 % from 550 to 720 nm. The 
pixel size is 9 jiva and the total CCD size is 3.7 cm. The 
pixel response non-uniformity is around 0.5 %. Pixels are 
coded on 16 bits, implying a dynamic range of 65535 ADU. 
The gain is 2 e-/ADU. A filter whose transmission starts 
at 600 nm is placed before the camera to eliminate blue 
light. Given the CCD quantum efficiency, the overall trans- 
mission (600 to 900 nm) is equivalent to that of a large R 
band. We use a CM 8 equatorial mount from Losmandy. A 
thermalized enclosure is used to avoid temperature fluctu- 
ations. The sides of this enclosure are made with wood and 
polystyrene. A double glass window reduces temperature 
variations and its accompanying turbulence on the optical 
path. Windows are fixed together by a teflon part and sep- 
arated by a 3 mm space filled with nitrogen to avoid vapour 
mist. The enclosure is thermalized to — 20°C and fans are 



used for air circulation. The ASTEP South instrument is 
shown at Dome C in figure [TJ 

In order to characterize the quality of Dome C for pho- 
tometric observations, we have to avoid as much as possible 
instrumental noises and in particular jitter noise, leading to 
a new observation strategy: the instrument is completely 
fixed and points towards the celestial South pole contin- 
uously. This allows also a low and constant airmass. The 
observed field of view is 3.88 x 3.88°^, leading to a pixel 
size of 3.41 arcsec on the sky. This field contains around 
8000 stars up to magnitude R = 15. This observation setup 
leads to stars moving on the CCD from frame to frame and 
to a widening of the PSF (Point Spread Function) in one 
direction, depending on the exposure time. 

Test observations were made at the Calern site 
(Observatoire de la Cote d'Azur) observing the celestial 
North pole, in order to choose the exposure time and the 
PSF size. A 30 second exposure time and a 2 pixel PSF 
FWHM (Full Width Half Maximum) lead to only 2 satu- 
rated stars and a limit magnitude around 14 (from Dome C 
the limit magnitude is increased to 15). An analysis of 
the celestial South pole field from the Guide Star Catalog 
GSC2.2 with these parameters taking into account the ro- 
tation of the star during each exposure leads to less than 
10 % of blended stars. Therefore we adopted these param- 
eters. 

Software programs were developed by our team to con- 
trol the camera, to run the acquisitions and to transfer and 
save the data. The instrument was set up at the Concordia 
base in January 2008. 




Fig. 1. ASTEP South at Dome C, Antarctica, January 
2008. 



2.2. The South pole field 

The distribution of stars in our field of view is shown in fig- 
ure [2l From the GSC2.2 catalog, we find nearly 8000 stars 
up to our limit magnitude of 15. We also simulate stellar 
populations in a field of 3.88 X 3. 88°^ centered on the celes- 
tial South pole us ing the Besancon model of the Galax}{3 
(Ro bin et alll2003D for R-band magnitudes between 10 and 
18 to calculate the dwarf ratio in the field. The comparison 
shows that the Besancon model overestimates the number 



^ ,http : //bison . obs-besancon . f r/modele/ | 
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Table 1. Number of stars in the 3.^ 
pole field. 
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of stars in the field by a factor ~ 2. However, we believe 
that the ratio of dwarfs to the total number of stars is, by 
construction of the model, better estimated. The bottom 
panel of fig. [2] shows that most of the stars brighter than 
magnitude R = 12 are giants (or more accurately larger 
than twice our Sun). 

Table [1] details the number of stars per magnitude range; 
the total number of stars is obtained from the GSC2.2 cata- 
log and the number of dwarfs is estimated using the relative 
fractions from the Besancon model. From magnitude 10 to 
15 we have 73.6 % of dwarf stars with radius R < 2Rq. This 
ratio is higher than in other typical fields used in the search 
for transiting planets such as C arina. Based on CoRoTlux 
simulations (Fressin et aI][2QQ7[ ). we expect that about one 
F, G, K dwarf in 1100 to 1600 should harbor a transiting 
giant exoplanet. The South pole field observed by ASTEP 
South is thus, in principle, populate d enough for the de - 
tection of transiting planets (see also lCrouzet et al1l2009l ). 
We will come back to a realistic estimate of the number of 
detectable exoplanets in section [5l 

The advantages of the South pole field are hence of 
course a continuous air mass, a high ratio of dwarfs to giant 
stars and a very low contamination by background stars. On 
the other hand, the field is less dense than regions closer to 
the galactic plane, so that the actual number of transiting 
planets in the field is smaller. 




14 

Magnitude 

Fig. 2. Top panel: Cumulative distribution of the number 
of stars in the South pole field as a function of their magni- 
tude in the R band. The plain line shows results from the 
Besancon model. The dashed line indicates output from the 
GSC2.2 catalog. Bottom panel: Ratio of dwarf stars with 
selected radii (less than 2, 1 and 0.5 i?©, respectively, as 
labeled) to the total number of stars in the South pole field 
as a function of R magnitude. 



2.3. Temperature conditions 

The instrument was set up during the Dome C 2008 summer 
campaign. The external temperature varied at this time 
between —20 and — 30°C. It was let outside without thermal 
control until the observations started at the end of April. 
In winter the external temperature varies between —50 and 
— 80°C. During the observations, the thermalized box is set 
to a temperature of — 20°C and the CCD to — 35°C. Because 
of self-heating, the electronics of the camera is around +5°C 
with some variations (see 13. 7p . 

3. Preliminary data analysis 

ASTEP South generates around 60 gigabytes of data per 
day. Since internet facilities at Dome C are limited to a 
low stream connection only few hours a day, a whole data 
transfer is impossible. Data are stored in external hard disks 
and repatriated at the end of the winter-over, leading to at 
least a 6 month delay between the observations and a full 
data analysis. We thus developed a software program for 
on-site preliminary data analysis, in order to have a day-to- 
day feedback of the observations. We detail here the results 
of this preliminary analysis. 



3.1. Preliminary data analysis software program 

We developed a software program running on the data at 
Concordia. For each image of a given day the mean inten- 
sity is computed. We then process only the 1000 x 1000 
pixel central part of the frame (0.95 x 0.95°^) for faster cal- 
culations. First, a point source identifier gives the number 
of detected stars and their location on the CCD. The 200 
brightest stars are matched to the CSC 2.2 catalog using a 
home-made algorithm, in order to identify the South pole 
on the CCD. The 30 brightest stars are fitted with a gaus- 
sian to derive the PSF size. Last, basic aperture photometry 
is performed for a set of 10 stars without any image calibra- 
tion. The identification of point sources, the gaussian fit and 
the aperture p hotometry use an IDL version of DAOPHOT 
(IStetsonlll987h . A point source is considered as a star if its 
fiux is 5 times larger then the sky noise. Aperture photom- 
etry is made with large apertures of diameter 12 and 20 
pixels, allowing to get all the fiux for bright stars. Although 
these large aperture are not adapted to faint stars, the low 
crowding in our field allows to get reasonable light curves. 
Of course this will be optimized during the complete anal- 
ysis of data. The camera and CCD temperature are also 
recorded. A small size binary file with these results is sent 
everyday by email. Plots shown in the following are in UTC 
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time as recorded by the software program (local time at 
Dome C is UTC +8). 



3.2. Magnitude calibration 

In order to convert ADU into magnitudes, we perform a 
preliminary magnitude calibration: we measure the flux of 
the stars on a typical image taken under dark sky and con- 
vert them into instrumental magnitudes. We then compare 
these magnitudes to the ones from the GSC2.2 catalog and 
obtain the so-called zero point. The image used is a raw im- 
age, but the local background including bias is subtracted 
when calculating the flux of each stars. 

We estimate that the error on these magnitudes should 
be ±0.3 magnitudes or less. First a comparison of the result 
for all the stars in a given image to that obtained with only 
the stars in the 1000 x 1000 pixel central part yields a 0.2 
magnitude difference. We estimate that the absence of a 
flat-field procedure is responsible for that difference and 
that its impact on our inferred sky brightness magnitude 
should be smaller. Second, while one may estimate that the 
GSC2.2 errors on the magnitudes of individual stars can be 
as large as 0.5, the large number of stars (~ 7000) implies 
that the mean error should be quite smaller. A 0.3 error on 
the inferred magnitudes hence appears to be a conservative 
estimate. 

In what follows, we will use this ADU to magnitude con- 
version only for the noon and full- moon sky brightness, not 
for the dark sky. This is because our preliminary analysis 
is based on data processed on the fly in Concordia which 
have not been de-biased. Variations in the bias level are of 
the order of 40 ADU. Given that uncertainty, we estimate 
that any measurement of magnitudes larger than 18 may 
have a bias error larger than 0.3 magnitudes and therefore 
refrain from mentioning those. 

A refined analysis of the full ASTEP South data with 
all available data is under way and will include an accurate 
de-biasing and magnitude calibration. 



3.3. Influence of the Sun 

We first consider the influence of the Sun on the photome- 
try. It is important to notice that although the Sun disap- 
pears below the horizon from May 4 to August 9, the sky 
background is always higher each day in the period around 
noon which is therefore less favorable for accurate photo- 
metric measurements. The minimum altitude of the Sun at 
noon occurs on June 21 and is 8.5° below the horizon. The 
height and width of the peak of intensity are the smallest 
around the winter solstice and increase before and after this 
date (figure [3]) . The increase is not linear but varies from 
one day to another, as also observed with the sIRAIT in- 
strument (^Strassmeier et al. 200S). We attempted to check 
whether this may be due to high altitude clouds but no 
correlation was found between the sky brightness and the 
quality of the night derived by studying the number of de- 
tected stars (see section HJ. 

Figure [4] shows variations of the mean intensity as a 
function of time for 3 clear days: June 22, July 20 and 
August 20. On June, 21 the height is typically 1600 ADU 
and images are affected during 4 to 6 hours. From our cal- 
ibration this corresponds to a magnitude of 15.3arcsec~^ 
in the standard R band. The residual noise calculated from 



the actual number of photons received from the sky in an 
aperture of 20 pixels (corresponding to a radius equal to a 
FWHM of 2.5 pixels) is 4 x 10~^. For larger apertures the 
noise will be smaller. Therefore this effect will have a mod- 
erate impact on the photometry. In July the height grows 
to typically 20000 ADU, i.e. magnitude 12.6 arcsec"^, and 
a noticeable brightness increase lasts for 7 to 9 hours. In 
August, this brightness increase lasts 9 to 12 hours. 




10 13 30 25 
Augygt 



10 i& 30 25 30 
September 



Fig. 3. Image intensity for all the winter. The bias level is 
around 2480 ADU and the scale is half of the CCD dynam- 
ics. Peaks around noon will affect the photometry. 
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5 10 15 20 

Time from OhOO (hours) 

Fig. 4. Image intensity during 3 typical days in June, July 
and August. The bias level is around 2480 ADU and the 
scale is half of the CCD dynamics. Due to the Sun, the 
sky background increases around noon. The corresponding 
magnitude for June 22 is 15.3 mag / arcsec^. 



The mean intensity of each image and the number of 
detected stars are plotted against the height of the Sun 
in figure [5l The fact that the sky intensity drops to an 
undetectable level when the Sun is below —13 ° appe ars to 
be in line with the result from 'Moore e t al.l (1200 8*) that 
the Dome C sky may be darker than other sites. However, 
this conclusion is at most tentative due to the absence of 
a bias removal and dark sky magnitude determination. We 
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notice that the R-band sky magnitude averaged over all 
observations for a Sun altitude of —9° is 16.6 arcsec"^, very 
similar to that obtained close to the zenit h for Paranal in 
the R-band, i.e. 16 to ITarcsec" -2 (iPatat et all l2QQ6l . see 
their figure 6). 



A typical increase of 80 ADU during the full Moon leads 
to a sky brightness of ~ 18.1 mag/ arcsec^. As discussed 
in section 13. 2[ this magnitude estimate may change by a 
fraction of a magnitude with a precise bias subtraction. 




5 10 15 20 25 30 




-12 -10 -S 
Sun oltitude (degree) 



Fig. 5. Mean intensity of each image and number of stars 
as a function of the Sun altitude when it is between 
-15° and -5°. The bias level is around 2480 ADU and the 
intensity scale is the whole CCD dynamics. Points: individ- 
ual measurements, blue and red: average during ascending 
and descending periods respectively. At 9° below the hori- 
zon the sky brightness is 16.6 mag / arcsec^. 



3.4. Influence of the Moon 

The influence of the Moon is shown in figure [6l The Moon 
is full on June 18, July 18, August 16 and September 15. 
An increased sky background is clearly seen around these 
days, up to 80 ADU in June, 100 ADU in July, 500 ADU in 
August and 70 ADU in September. The full Moon in June 
and July corresponds to a good weather, without clouds, 
and the increase in intensity is low enough to allow pho- 
tometric observations. In contrast, during the full Moon 
of August the weather was very bad with high temperature 
(up to — 30°C), strong wind at ground level (up to 11 m/s), 
and a very cloudy sky. The very high background is thus 
interpreted as due to the reflection of moonlight by clouds. 



Fig. 6. Image intensity. The bias level is around 2480 ADU 
and the scale is 1 % of the CCD dynamics. The sky back- 
ground increases during 10 days around the full Moon up 
to typically 18.1 mag/ arcsec^. 



3.5. Point Spread Function variations 

PSF variations are a crucial issue for photometry. We inves- 
tigate here the PSF variations in the ASTEP South data. 
For each image, the 30 brightest stars are fitted with a gaus- 
sian PSF and their FWHM in both direction is calculated 
using DAOPHOT. The mean of the FWHM across the en- 
tire image is shown as a function of time in figure [71 This 
mean FWHM varies between 1.5 and 3.5 pixels over the 
winter. 

Two kinds of variations are present. First, PSF vari- 
ations on a timescale smaller than one day are ob- 
served. We compare them to independent seeing measure- 
ments at Dome C provided by three dedicated Differential 
Image Motion Monitors (DIMM), two of them forming a 
Generalized Seeing M onitor (GSM) (fo r a description of 
these instruments see IZiad et al.|[2008D . In order to con- 
sider only the PSF variations of period smaller than one day 
we subtract to the FWHM the difference between the me- 
dian FWHM and the median seeing for each day. Figure [8] 
shows that on this day timescale the corrected FWHM 
and the seeing are clearly correlated: the PSF variations 
on short timescales are mostly due to seeing variations. As 
previously discussed, the seeing at the ground level where 
ASTEP South is placed is rather po or: the median value in 
winter at 3 meters high reported bv lAristidi et al.l (|2009l ) is 
2.37 arcsec with stability periods of 10 to 30 minutes. This 
explains the short-term variations of our PSFs. 

On a timescale larger than one day the correlation is not 
true anymore. This shows that another cause of PSF varia- 
tions is present. For this larger timescale, two regimes seem 
to be present, one with a PSF around 1.5 pixels and another 
with a PSF around 3.0 pixels. These variations are associ- 
ated with an asymmetrical deformation of the PSF. This 
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suggests an instrumental cause of PSF variations such as 
astigmatism and decolHmation. Indeed, temperature inho- 
mogeneities in the thermal enclosure cause mechanical and 
optical deformations. Unfortunately these large timescale 
variations prevent us from estimating the seeing at the 
ground level directly from our photometric data. 
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Fig. 7. Mean values of the size of the stars (FWHM) on the 
CCD in pixels (top panels) and their asymmetry (bottom 
panels) as a function of the observing day for the months of 
June (top left), July (top right), August (bottom left) and 
September (bottom right). In the top panels, the blue and 
green curves correspond to the values of the FWHM in the 
X and y directions, respectively. The mean FWHM values 
are obtained through a spatial average on the CCD. 
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Fig. 8. Correlation between the PSF FWHM and inde- 
pendent seeing measurements at Dome C, for timescales 
smaller than one day (the FWHM is corrected for larger 
timescale variations). Direct seeing measurements from 
three Differential Image Motion Monitors are shown in blue, 
green and red. A linear regression gives a slope a = 0.59 
with a correlation coefficient r = 0.65. The PSFs are clearly 
affected by seeing variations at the ground level. 



3.6. Astrometry and pointing variations 

Ideally the pointing should remain stable during all the 
winter, meaning that the South pole must stay at the same 
place on the CCD. The position of the South pole on the 
CCD is found on each image using a home-made field- 
matching algorithm. The precision of this algorithm is typ- 
ically 0.2 pixels. The results for a typical day and for all 
the winter are shown in figure [9) 

First we have a variation of this position with a period of 
one sidereal day. This is due to an incomplete correction of 
astrometric effects. Indeed, the star coordinates from the 
GSC2.2 catalog were corrected only for the precession of 
the equinox from the J2000 epoch to January 1, 2008. The 
remaining error on the star coordinates led to an error of 
10 pixels (34 arcsec) in the determination of the position of 
the pole. We then corrected the GSC2.2 coordinates from 
the precession of the equinox using the real observation 
date, and from the nutation and the aberration of light (or 
Bradley effect). After these astrometric corrections the pole 
stays within 2 or 3 pixels during the day. 

Second, the pole is drifting during the winter. The am- 
plitude is 10 pixels (34 arcsec) in 50 days, from June 12 to 
July 31. From the orientation of the CCD we find that this 
drift is oriented vertically towards the North. This may be 
due to mechanical deformations of the instrument, atmo- 
spheric changes or a motion of the ice under the instrument. 
In any case this effect is very small. 



3.7. Camera temperature variations 

The CCD is cooled down to — 35°C without any variations. 
In contrast the electronic part of the camera oscillates be- 
tween +4 and +8°C with a one hour period (figure [TQ|) . 
A threshold effect explains these variations: the thermal- 
ized enclosure is not heated continuously but only when it 
passes below a threshold temperature. A direct consequence 
is seen on the bias images. The bias level oscillates with the 
same period and an amplitude of 10 ADU. The mean inten- 
sity of science images is affected in the same way. The bias 
level is plotted against the camera temperature in figure [TTl 
and shows a hysteresis behavior. For a given temperature, 
the bias level is lower if the temperature is increasing than 
if the temperature is decreasing. The hysteresis amplitude 
is around 5 ADU. An explanation can be that the tempera- 
ture sensor is not exactly on the electronics but is stuck on 
a camera wall which may be sensitive to temperature vari- 
ations with a time lag compared to the electronics. It may 
also be due to the electronics and sensor having different 
thermal inertia. 



4. Duty cycle 

A main objective of ASTEP South is to qualify the duty 
cycle for winter observations at Dome C. The observation 
calendar for the whole 2008 campaign is shown in figure [T2l 
April and May were mainly devoted to setting up the in- 
strument and software programs. Continuous observations 
started around mid- June. Since then, very few interruptions 
occurred and data were acquired until October. The effect 
of the Sun and of the Moon has already been discussed 
in section [3l We present here some technical limitations to 
the duty cycle, and quantify the photometric quality of the 
Dome C site for this campaign. 



Crouzet et al.: ASTEP South 



7 



2035 



2030 



^ 2025 



2020 



2015 




2010 2015 2020 2025 2030 2035 
X (pixel units) 



2035 



^ 2030 
cn 

*E 

2025 
2020 



2015 



1 1 1 1 1 1 1 1 1 1 1 1 

~ 24/06 02/07 
■ 1^ ++ 12/06^ 


1 1 1 1 1 1 1 1 1 1 1 1 1 

T 31/07 


w 

. 1 .... 1 .... 1 . 


* s- 

11/07 ■ 
... 1 .... 1 ... . 



2010 2015 2020 2025 
X (pixel units) 



2030 2035 



Fig. 9. Position of the pole on the CCD. Top: 9 images on 
July 11 after various astrometric corrections (blue: approx- 
imate correction of the precession of the equinox, green: im- 
proved correction of the precession of the equinox, brown: 
correction of the nutation, red: correction of the Bradley ef- 
fect). Bottom: 7 days between June 12 and July 30, showing 
a drift over the winter (dark blue: June 12, light blue: June 
24, green: July 2, yellow: July 11, orange: July 17, red: July 
30). 



ments should be equipped with converters to avoid pos- 
sible damages. 

— The instrument is submitted to temperature gradients 
inside the thermal enclosure, and to the external tem- 
perature during power cuts or when opening the box. 
This leads to mechanical constraints resulting in decol- 
limation and astigmatism. 

4.2. Weather conditions at Dome C 

A first experiment to measure the winter c l ear sky frac- 
tion at Dome C was made by lAshlev et al.l (|2005at ) with 
ICECAM, a CCD camera with a lens of 30° field of view. 
Every 2 hours from February to November 2001, ten im- 
ages of the sky were taken and averaged. An analysis of all 
the images yielded a fraction of 74 % of cloud- free time. An 
analysis by iMosser fc Aristidil (|2QQ7[ ) for the 2006 winter 
yielded an estimate of 92 % of clear sky fraction by re- 
porting se veral times a day th e presence of clouds with the 
naked eye. iMoore et al.l (|2QQ8l ) derive a clear sky fraction of 
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Fig. 10. Camera temperature, bias level and image inten- 
sity on July 11. The camera temperature varies between +4 
and +8°C with a period of one hour and affects the bias 
level and the image intensity to about 10 ADU. 



4.1. Technical issues 

Technical issues encountered during this campaign limited 
the duty cycle. We show here typical issues that instru- 
ments at Dome C have to face with. We believe these can 
be mostly overcome with appropriate technical solutions. 

— First, the shutter did not close and got damaged at tem- 
peratures below 0°C. We had to change the shutter and 
build a special thermalization device to warm it. 

— In order to install the camera again after changing the 
shutter, the thermalized box was opened and suddenly 
cooled by the ambient air at ~ — 60°C. As a result, 
cables not made in teflon broke as well as the camera 
USB connection. These had to be replaced. 

— Outside instruments are affected by power cuts lasting 
for a few minutes to a few hours. The fraction of time 
lost for observations is negligible, however next instru- 
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Fig. 11. Bias level against the camera temperature for July 
11. The blue points correspond to an increasing tempera- 
ture and the red points to a decreasing temperature, show- 
ing a hysteresis behavior. 
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Fig. 12. ASTEP South observation calendar for the 2008 
campaign. April and may were mainly dedicated to solving 
technical problems. Continuous observations started mid- 
June with very few interruptions until the end of the winter. 



79 % for the winter 2006 from the Gattini instrument using 
the number of stars and the extinction across the images. In 
a previous work, we derived a clear sky fraction of 74 % for 
the 2008 winter from the ASTEP South data, considering 
that the sky is clear if we h ave more than half o f the stars 
detected on the best images (jCrouzet et al.|[2QTQl ). However 
this result is dependent on this ad hoc criterion. We reeval- 
uate here this fraction by avoiding such an arbitrary limit. 

4.2.1. Method 

A new measurement of the clear sky fraction is made with 
ASTEP South using a method sensitive to thin clouds, 
based on the number of stars detected in the field. In or- 
der to do so, we need to evaluate the number of stars that 
should be detected on any given image if the weather was 
excellent. Our PSF size varies due to fluctuations of the see- 
ing and of the instrument itself, and the background level 
also changes due to the presence of the Moon and the Sun. 
Since these are not directly related to weather, we need 
to derive how the number of detectable stars changes as a 
function of these parameters. (Note that thin clouds should 
affect the seeing in some way, however a posteriori exami- 
nation of the data shows that this effect is small compared 
to the global attenuation due to clouds). 



4.2.2. Identifying point sources with DAOPHOT 

The 1000 X 1000 pixel sub-images contain up to 500 stars 
of varying magnitude. Our automatic procedure for find- 
ing point sources uses the FIND routine from DAOPHOT. 
Specifically, in this procedure a star is detected if the central 
height of the PSF is above the local background by a given 
number of standard deviations of that background. This 
threshold parameter a is chosen by the user. We choose 
a = b. 

4.2.3. A model to evaluate the expected number of stars 

The full width half maximum co of our PSFs is typically 
2.5 pixels. To evaluate whether a star is detected or not we 
compare the amplitude A of the PSF to the noise in the 
central pixel. We consider two kind of noises: the photon 
noise from the sky background N^ky = y^Fsky and the read- 
out noise A^ron- The noise in the central pixel is hence N = 

'\/^sky + ^ron • ordcr to obtaiu A as a function of co we 
consider a gaussian PSF. The amplitude A of a gaussian 
is A = cF/cj'^ with F beeing the total flux under the PSF 
and c ~ 0.88. The condition to detect a star A > aN can 
thus be expressed as F > aNuj'^/c. The limit magnitude m 
is therefore: 

m = -2.5 \og{aNuj'^/c) + Z (1) 

with Z being the photometric calibration constant. We have 
^ron = 10 electrons and set a = 5 as we do in DAOPHOT, 
and since the instrument is not calibrated photometrically 
we use as ad hoc constant Z = 21.6. As an example, typical 
values in our data are Fgky = 80 electrons and a; = 2.5 
pixels. This yields m = 14.9. 

To derive the number of stars A/"* expected in a 
1000 X 1000 pixel sub- image from this limit magnitude we 
use a typical image taken under excellent weather condi- 
tions. We calculate the distribution in magnitude of the 
detected stars and fit it with a 3^^ order polynomial. For 
magnitudes larger than 14 the number of stars increases 
more slowly because they are becoming too faint to be all 
detected. We therefore extend the fitting function with a 
constant slope. The following relation provides our assumed 
number of stars as a function of the limit magnitude: 



( asm^ -\- a2m'^ -\- aim -\- ao if m < 14 . . 
* " \log7V*i4 + 0.2 (m- 14) if m > 14 

where as = 0.013, a2 = —0.664, ai = 11.326, 
ao = —61.567 and A^*i4 is the number of stars detected for 
m — 14. Equations ([1]) and ([2]) thus provide the number of 
stars that should be detected for a clear sky given a value 
of sky background and FWHM. In order to test the validity 
of this relation, we compare this to the maximum number 
of stars detected in our images for given values of FWHM 
and sky background. (By choosing the maximum number 
of stars, or more precisely the number of stars which is ex- 
ceeded only 1 % of the time, we ensure that we consider 
only images taken under excellent weather conditions). We 
find that both agree with a standard error of 6.6 % and a 
maximum error of 15 %. 



Crouzet et aL: ASTEP South 



9 



4.2.4. Comparison to the measured number of stars 

We use this model to compare the measured to the ex- 
pected number of stars for each data point, given its FWHM 
and sky background. The resulting distribution in figure [13] 
shows two features: a main peak centered around 1 with 
most of the points, meaning that for these points the 
weather is excellent, and a tail for which the measured 
number of stars is smaller than expected, corresponding 
to veiled weather. 

The spread around the main peak of the histogram is 
partly due to measurement errors in the sky background 
and FWHM. The limit between this natural spread and 
the veiled weather data points must be defined to calcu- 
late the weather statistics. To do this we use a large set of 
sky background and FWHM values coherent with our mea- 
surements, and calculate the expected number of stars for 
each point (Fsky,cc;) according to our model. We then add 
some random errors to this set of values and calculate again 
the expected number of stars. For each point {Fsky^co) we 
compare the expected number of stars obtained with the 
added errors to the one without errors. This yields a theo- 
retical distribution of the number of stars that accounts for 
measurement errors. 

The sky background value is typically around 40 ADU. 
We estimate the error to be around 10 % of this value, i.e. 
only 4 ADU. This error is not dominant and the spread of 
the theoretical distribution is mainly due to the error on the 
FWHM. The FWHM is typically between 2 and 3 pixels, 
and we suppose again an error of 10 %, i.e. 0.25 pixels. 
However the shape of the theoretical distribution does not 
fit perfectly to the data. Instead we find that we can fit the 
two sides of the main peak in figure [13] with two different 
errors on the FWHM, corresponding to a low and a high 
measurement error, respectively 0.17 and 0.25 pixels. 

The data points fitting into the theoretical distribution 
with the low error are considered as excellent weather (red 
part in figure [T3|) . Those between the low and high error 
distributions can be either due to a large measurement error 
or to veiled weather, thus they are considered as uncertain 
weather (orange part). The data points outside the low er- 
ror distribution correspond to veiled weather (yellow part) . 

4.2.5. Weather statistics for the winter 2008 

This analysis gives a fraction of time between 65.2 % and 
79.2 % with excellent weather and between 20.8 % and 
34.8 % of veiled weather. Only the periods with data when 
at least few stars are visible are considered here, exclud- 
ing in particular the white-out periods. During the winter 
the acquisitions were stopped during 13.7 % of time be- 
cause of very bad weather, so the previous numbers must be 
multiplied by 1 — 0.137 = 0.863. The weather statistics for 
the 2008 winter at Dome C are therefore: between 56.3 % 
and 68.4 % of excellent weather, 17.9 to 30 % of veiled 
weather during which stars are still visible and 13.7 % of 
bad weather (figure [T4|). For comparison the fraction of 
photometric weather during night-time is 62 % at La Silla 
and 75 % at Paranal as provided by M.Sarazirj^ (see also 
lArdeberg et al.l 19901) , whereas it is only 45 % at Mauna 
Kea (Ortolani" "2003', and references therein) though more 
recent results report 56 % (Steinbring et al...2009.) . Dome C 
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Fig. 13. Ratio between the measured and the expected 
number of stars (black line). Theoretical histograms tak- 
ing into account low and high measurement errors allow 
to identify 3 parts in the data: excellent weather in red, 
uncertain weather in orange and veiled weather in yellow. 



is therefore very competitive compared to other observing 
sites. 
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Fig. 14. Cumulative histogram of the ratio between the 
measured and the expected number of stars. The data 
are distributed between excellent, uncertain and veiled 
weather. The white part is the fraction of time during which 
the acquisitions where stopped because of very bad weather 
(white-out). 



4.2.6. Validation 

To validate the method, we compare the number of stars 
to the intensity received from the stars, as both should be 
smaller if clouds are present. We measure the intensity re- 
ceived from nine stars of magnitude 8 to 12 and sum them 
after dividing each star by its median value. We then nor- 
malize this sum by the maximum value, more precisely by 



|http : //www . eso . org/gen-f ac/pubs/astclim/paranal/cloudd yteaitteanurfll.^^ 1 % highest values. We use only the periods 
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with a moderate sky background i.e. when the Sun is be- 
low — 13° and excluding the full Moon periods (the result is 
however very similar using all data points). Figure [15] shows 
the normalized star intensity as a function of the ratio be- 
tween the measured and the expected number of stars. As 
expected, both parameters are directly correlated, thus val- 
idating the method. We further note that for data points 
taken under excellent or intermediate weather conditions, 
i.e. when N^/N^th > 0.7, the normalized intensity received 
from the stars is 0.77 ±0.22. Both parameters are thus good 
indicators of the cloud cover in the field of view. 
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Fig. 15. Normalized intensity of a selected sample of nine 
stars as a function of the ratio between the measured and 
the expected number of stars based on theoretical estimates 
(see text). 



4.3. Auroras 

Aurorae were first feared to be a source of contamination 
for long-term photometric data. However, it is to be noted 
that Concordia is a favorable site in that respect: auroras 
occur mainly in the auroral zone, a ring-shaped region with 
a radius of approximately 2500 km around the magnetic 
pole. The Concordia site is located well within this ring, 
only 1300 km from the South magnetic pole, and thus has a 
much lower auroral activity than other sites (e.g. Dome A). 

During the winter-over 2008, the Concordia staff re- 
ported 2 auroras on July 30 and 31. On July 30, a careful 
examination of the ASTEP South data indicates a possibil- 
ity of auroral contamination between 14:12 and 14:24 UTC, 
but it cannot be distinguished from thin clouds. The July 
31 data were unfortunately lost during the copy, probably 
because of a hard drive glitch (the only instance of that oc- 
curring) so that we could not attempt to check the images 
for that day. 

In any case, the ASTEP South 2008 data were not 
contaminated by auroras, confirming the low contribu- 
tion of aurora s to t he sky brightness as suggested by 
iDempsev et al.l (j2005f ). It will be interesting to see whether 
it remains true when progressively moving towards a max- 
imum solar activity in 2012. 



4.4. Observing time and photometric quality of Dome C 

The duty cycle for the 2008 campaign of ASTEP South 
is represented in figure [T6l The limit due to the Sun, the 
observing time and the excellent and intermediate weather 
fractions are shown for each day, as well as the white-out pe- 
riods. We acquired 1592 hours of data with ASTEP South 
on a single field during the 2008 campaign. From the pre- 
vious analysis we have 1034 hours with excellent or uncer- 
tain weather. As a comparison, simula tions based on the 
method described in lRauer et al.l (|2008l ) show that the time 
usable for photometry in one year at La Silla for the field 
with the best observability is typically 820 hours per year 
(see section [5] for more details). Moreover, the white-out 
periods at Dome C last typically from one to a few days, 
allowing extended periods of continuous observations be- 
tween them. For example we observed every day during 
one month between July 9 and August 8. Considering the 
excellent and uncertain weather and the hours lost because 
of the Sun, the fraction of time usable for photometry for 
this one month period is 52 %. In La Silla the fraction of 
time usable for photometry for all one month periods be- 
tween 1991 and 1999 has a mean of 27 % with a maximum 
of 45 % in April 1997 (from the La Silla weather statis- 
tic^ multiplied by the night-time fraction). This shows the 
very high potential of Dome C for continuous observations 
during the Antarctic winter. 
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Fig. 16. Daily observing time fraction for ASTEP South 
in 2008 as a function of the observation period. The light 
blue and dark blue regions indicate the fraction of time for 
which the Sun is lower than —9° and —13° below the hori- 
zon, respectively. Periods of excellent, uncertain and veiled 
weather as observed by ASTEP South are indicated in red, 
orange and yellow, respectively. White areas correspond to 
periods during which observations were not possible, either 
because of the Sun altitude or because of bad weather. 



5. Planet detection probability 

As shown by iPont fc Bouchvl (|2005l ) , the high phase cover- 
age at Dome C should improve the efficiency of a transit 
survey. Here we investigate the potential of ASTEP South 



|http : //www . Is . eso . org/sci/f acilities/lasilla/astclim/weather/l 
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for transit detection using CoRoTlux and compare the 
Dome C and La Silla observing sites. CoRoTlux performs 
statistical simulations of transit events for a survey given 
the star distribution in the field of view, t he instrumen- 
tal p arameters and the observation windows (|Fressin et alJ 
l2QQ7l [2009). In all simulations, the star distribution is the 
one of the South pole field (see section [2^ . We use the 
GSC2.2 catalog for stars from magnitude 10 to 14.5 com- 
pleted with a distribution from the Besancon model up to 
magnitude 18 (after scaling it to the GSC2.2 catalog for 
low magnitude stars). The target stars range from magni- 
tude 10 to 15 and the background stars from magnitude 
15 to 18. The instrumental parameters are always those of 
ASTEP South. We perform three simulations correspond- 
ing to three survey configurations differing only by their 
observation windows. 

The first set of observation windows that we have 
used in our simulations corresponds to the periods during 
which ASTEP South actually ran in excellent or uncertain 
weather conditions in 2008, i.e. the red and orange parts 
of the duty cycle in figure [161 This provides us with the 
potential yield of the 2008 campaign in terms of detections 
of transiting planets. 

We also want to compare Dome C and La Silla. 
In that purpose we consider an ideal campaign for an 
ASTEP South-like instrument for which the observation 
windows are determined only from the altitude of the Sun 
and weather statistics at that site. For Dome C, we apply 
the weather statistics presented in section l4.2.5l to an entire 
winter season in order to generate the second set of obser- 
vation windows. For this second simulation, we incorporate 
over the Sun limited duty cycle 13.7 % of white-out periods 
and 17.9 % of randomly distributed cloudy periods lasting 
less than one day. 

For La Silla, we generate a third set of observation win- 
dows using the monthly weather statistics acquired from 
1987 to 2007 [3. The weather statistics for each month is 
taken as the mean of the photometric fraction for this given 
month over all years. At La Silla, one cannot simply stare 
at the South pole field continuously because it is low over 
the horizon. A best pointing can be found that maximizes 
the observation time as a combination of weather statis- 
tics, night-time an d airmass (for a complete description 
of the method see iRauer et al.|[2008[ ). The resulting field 
with the best observability is centered on RA = 18h30' and 
DE= — 58°54'. For consistency with the other simulations 
we use the same stellar population as for the South pole 
field, and consider that photometric observations are pos- 
sible when the Sun is less than —9° below the horizon. The 
resulting duty cycles for weather and Sun limited observa- 
tions of a single field over one year are shown in figure [T71 
for both Dome C and La Silla. The total observing time is 
typically 2240 hours for Dome C and 820 hours for La Silla. 

A large number of runs (~ 3000) are performed for each 
survey configuration in order to have a significant statistic. 
The results of the simulations provide the number of de- 
tectable planets. We assume that only transiting planets 
with a signal to noise ratio higher than 10 are detectable. 
This yields 1.08 planets for the ASTEP South 2008 cam- 
paign (i.e. 3244 planets over 3000 runs), 1.62 planets for a 
whole winter at Dome C, and 1.04 planet for La Silla. These 
numbers are low because ASTEP South is a small instru- 



ment, however the number of planets is notably higher for 
a survey from Dome C than from La Silla. The resulting 
planet detection efficiency is shown in figure [181 The detec- 
tion efficiency is defined as the number of detectable planets 
divided by the total number of simulated planets. In spite 
of technical problems at the beginning of the winter, the 
detection efficiency for the ASTEP South 2008 campaign 
is equivalent to the one obtained for one year at La Silla. 
When comparing an instrument that would run for the en- 
tire observing season, the detection efficiency is found to 
be significantly higher at Dome C than at La Silla both in 
terms of planet orbital period and transit depth. For exam- 
ple we have an efficiency of 69 % at Dome C vs 45 % at La 
Silla for a 2-day period giant planet, and 76 % at Dome C 
vs 45 % at La Silla for a 2 % transit depth. The detection 
efficiency decreases for planets with longer orbital periods, 
but is even more favorable to Dome C relatively to La Silla. 
On the other hand, it is true that a mid-latitude site offers 
more available targets. However, we believe that this shows 
the high potential of Dome C for future planet discoveries. 
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Fig. 17. Observing time fraction as a function of observ- 
ing period at Concordia and La Silla. The blue and grey 
envelopes indicate values obtained by imposing a Sun al- 
titude lower than —9° below the horizon, for Concordia 
and La Silla, respectively. The red histogram is an exam- 
ple of a generated window function for Concordia using the 
weather statistics obtained from ASTEP South. The green 
histogram is generated for the field which is observable the 
longest with a high-enough airmass at La Silla and using 
the 1987-2007 weather statistics of the site. 



6. Conclusion 

ASTEP South, the first phase of the ASTEP project ob- 
served 1592 hours of data. Night-time photometric obser- 
vations started in a nearly continuous way around mid- 
June and proceeded to the end of September, when the sky 
became too bright even at midnight local-time. Our pre- 
liminary analysis showed that the Sun affects our photo- 
metric measurements when it is at an altitude higher than 
— 13° below the horizon. The sky brightness at dusk and 
dawn appears to vary quite significantly from one day to 
the other, but its mean is very similar to results obtained 



Ihttp : //www . eso . org/gen-f ac/pubs/astclim/paranal/cloud3(?k)Bettxidih£L Izfeaaith at Paranal (an R-band sky-magnitude 
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Fig. 18. Calculated efficiency of detection of transiting gi- 
ant planets for a single field as observed by an ASTEP- 
South like survey during a full season, and as a function 
of the orbital period (top) and transiting depth (bottom). 
Dotted blue lines: Detection efficiency for the ASTEP South 
2008 campaign. Plain blue line: Detection efficiency for a 
full winter at Dome C for a circumpolar field limited only 
by the weather statistics and the constraint of a Sun alti- 
tude below —9°. Plain red line: Same as before, but for a 
survey at La Silla and the field with the best observability 
over the year (see Fig. 17 and text). 

R = 16.6 arcsec"^ for a Sun altitude of —9°). The full Moon 
yields a sky brightness of R ~ 18.1 arcsec"^. Apart from one 
possible instance lasting only 12 minutes, auroras had no 
noticeable impact on the data. 

An identification of the stars in the field allowed us to 
retrieve the precise location of the celestial South pole on 
the images and show that the pointing direction is stable 
within 10 arcseconds on a daily timescale for a drift of only 
34 arcseconds in 50 days. On the basis of the number of 
identified stars and of a model to account for PSF varia- 
tions and sky brightness, we retrieved the weather statistics 
for the 2008 winter: between 56.3 % and 68.4 % of excel- 
lent weather, 17.9 % to 30 % of veiled weather (when the 
probable presence of thin clouds implies a lower number of 
detected stars) and 13.7 % of bad weather. 

An analysis of the yield of transit surveys with our 
weather statistics at Dome C compared to those at La Silla 



showed that the efficiency to detect transiting planets in 
one given field is significantly higher at Dome C (69 % 
vs. 45 % for 2-day period giant planets with an ASTEP 
South-like instrument in one season). The prospects for the 
detection and characterization of exoplanets from Dome C 
are therefore very good. Future work will be focused on a 
detailed analysis of the full ASTEP South images. The sec- 
ond phase of the project includes the installation of ASTEP 
400, a dedicated automated 40-cm telescope at Concordia 
and its operation in 2010. 
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